Collective excitations and instability of an optical lattice due to unbalanced pumping 
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We solve self-consistently the coupled equations of motion for trapped particles and the field of a one- 
dimensional optical lattice. Optomechanical coupling creates long-range interaction between the particles, 
whose nature depends crucially on the relative power of the pump beams. For asymmetric pumping, travel- 
ing density wave-like collective oscillations arise in the lattice, even in the overdamped limit. Increasing the 
lattice size or pump asymmetry these waves can destabilize the lattice. 
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Optical lattices (OL) are perfectly periodic arrays of parti- 
cles trapped by the standing wave interference pattern of sev- 
eral laser beams. They have important applications as model 
systems for solid state physics as well as for quantum infor- 
mation science. The back-action of the trapped particles on 
the trap light is carefully avoided in most OL experiments. 
However, it is known to give rise to intriguing phenomena in 
related systems, e.g., cavity cooling |Q]], mirror cooling J2l, 
and optical binding fH. For OL's this back-action has been 
predicted 10] and observed to reduce the lattice constant 
compared to the naive expectation. 

In this Letter we consider the dynamical effects of optical 
back-action in a one-dimensional OL, brought about by tun- 
ing a hitherto neglected parameter, the asymmetry in the in- 
tensities of the lattice beams. Due to the back-action the trap 
light mediates an interaction between the particles, which is 
substantially altered by this asymmetry. Net energy and mo- 
mentum flow is induced through the OL, relating it to crystals 
driven far from equilibrium, e.g., arrays of vortices in a type-II 
superconductor [7], and trains of water drops dragged by oil 
(H. The phonon-like traveling density waves characteristic of 
these systems become the elementary excitations of the OL as 
well, and can destabilize it, even in the presence of arbitrarily 
strong viscous damping. The excitations arise resonantly at 
specific values of the asymmetry, which allows for tuning the 
dispersion relation of the lattice. Moreover, the light-mediated 
interaction in the OL is of infinite range, and thus all these ef- 
fects depend heavily on the size of the lattice. As absorption 
inevitably leads to pumping asymmetry, this dynamic insta- 
bility limits the size of any near-resonant OL. 

We consider a dipole trap formed by two counter- 
propagating phase locked laser beams with frequency lo = ck. 
The waist of the trap is much larger than the wavelength 
A = 2ir/k, so the light field is essentially 1 dimensional along 
x. The two beams have unequal intensities: the electric field 
incident from the left is E{x) = Eoe tkx ~ luJt , from the right, 
E(x) = E x e- ikx - lut , with E x = e^VVEo, and the rela- 
tive phase. Besides the pump power ratio V > 1, we introduce 
another measure of the asymmetry: A = \Ei / E \ — \E / Ei\. 
We consider particles of linear polarizability a and mass 
rriA pre-cooled down to very low temperatures (possibly pre- 
trapped) and trapped by the dipole force in the light field. 
These can be atoms, the lasers being detuned to the red of 
a specific transition so far that it is not saturated and spon- 
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FIG. 1: (color online) A dipole trap created by two lasers of equal 
frequency but unequal power. The intensity (in light red), mirrored 
for better visibility, ranges between I m i n — i eo c ( | Eq \ — \ Ei | ) 2 and 
/ m „ = ieo c ([i?o + \Ei\) 2 . Trapped particles form disk-shaped 
clouds (in dark blue), and are modeled as beam splitters. Due to 
the pump asymmetry, the electric field has no nodes. Back-action of 
trapped particles distorts the field and reduces the lattice constant. 



taneous emission can be ignored. Alternatively, they can be 
plastic beads trapped in water, as in, e.g., JSJ 1 1011 . albeit of 
size well below A so that complications of Mie scattering are 
avoided. If the particles are cold enough, they gather at the 
antinodes, forming N disk-shaped clouds of axial size much 
smaller than A. For simplicity we assume that each cloud has 
the same number of particles, and thus identical surface den- 
sity ?y, surface mass density m — r/rriA, and dimensionless 
polarizability £ = krja/(2eo). The setup is sketched in Fig.Q] 
We now take the back-action of the particles on the light 
field into account. Following J3l, this is achieved by solving 
the scalar Helmholtz equation, with the N clouds represented 
by Dirac-(5 distributions of linearly polarizable material, 



(<9 2 + k 2 )E{x) = -2kE(x) V (5{x - Xj ) 
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Throughout this Letter we assume £ E M, neglecting spon- 
taneous emission and scattering into other transverse modes, 
justified as long as the laser beams are far detuned from any 
resonance of the trapped particles. Note that although these 
approximations can be relaxed by setting ( e C, very close to 
resonance the reabsorption of spontaneously emitted photons 
plays an important role in the dynamics ll ill , and this is not 
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easily incorporated into this model. 

The solution of Eq. (Q} between two clouds is a superposi- 
tion of plane waves, E(xj-\ < x < Xj) = Aje lk ( x ~ x ^ + 
B j e -ik(x-xj) _ (j.^ik^x-x,-!) + Dj^e-iKx-xj-i)^ The 

clouds constitute boundary conditions for the field: 

E(x = Xj -0) = E(x = Xj +0); (2a) 
d x E(x = Xj - 0) = d x E(x = Xj + 0) + 2kC,E{xj). (2b) 

This amounts to representing each cloud as a beam splitter 



(BS) at x 



with reflection and transmission coefficients 



r = </(l - <) and t = 1/(1 - <), so that C = -ir/t fl. 

Since E(x) is not differentiable at the cloud position, as 
given by Eq. d2bl and shown on Fig.Q] we need to calculate 
the dipole force on the cloud carefully. Integrating the force 
over a finite cloud and then taking the Dirac-5 limit, we obtain 



Fj = ^(d x E 2 (xj-0) + d x E 2 (xj 
o 



0)) 



(3) 



for the force on a unit surface of the cloud, averaged over an 
optical period. This formula can also be derived based on the 
amount of momentum transferred to the cloud by the field, via 
the Maxwell stress tensor, as in lfl2tl . 

For a single cloud at steady state, both Fj = and Eqs. (f2]i 
must hold, which is only possible if 



(A < 2. 



(4) 



This simple equilibrium criterion can be intuitively under- 
stood in the following way. If \E Q \ Z < \Ei\' ' , more photons 
are incident on the right of the BS than the left, giving a force 
on it. If enough light is transmitted (|t| > 1 \r\ A), and the in- 
terference is favourable (depending on the position of the BS), 
the imbalance in the outgoing number of photons is enough to 
counteract this force, leading to a steady state. 

For several clouds trapped by the same light, at steady state 
Fj has to vanish on all components of the system, which with 
Eqs. (f2]) and <[3j implies that E(x) and \d x E(x)\ are the same 
to the left and right of any component. As a result, |£? 2 (a;) | = 

\E \ 2 + \Ei | 2 + 2 \E Ei | cos(2fca; - everywhere in the 

sample, the clouds only contribute to the phase: $>(xj < x < 
Xj + i) — X^=i Xl< me phase slip at the /'th cloud depending 
on the polarizability Q of the cloud as 
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FIG. 2: (color online) The lattice constant as a function of the asym- 
metry is shown in thick (red) curves for £ = 0.01, f = 0.1, £ = 0.5, 
f = 1, £ = 2. Shaded (green) areas indicate regions of stability 
(see page 4), for N < 800 (darkest shade), N < 100, N < 10 and 
N < 2 (lightest shade). The white area is unstable, see Eq. {4}. The 
orange circle marks the parameter regime of Fig. [4] 



For .4 = this gives d symm = ^(1 — 2atan(£)/7r) as in 
|4j]. For a given £, increasing A causes the phase shift x to 
increase, and d to be reduced, as illustrated in FigfJ] (thick 
red lines). For A > 2/£, the inequality © is violated, the 
stronger beam pushes all the particles away. At A = 2/£ the 
lattice constant d is, remarkably, exactly half of d symm \ 



drain (C) = j- (t - 2 ata<) . 



(7) 



The fact that an equilibrium lattice configuration exists is 
only physically relevant if this equilibrium is dynamically sta- 
ble. The dynamics of an OL is given by 



mxj = —\iXj + Fj(xi..xisr), 



(8) 



cos xi (0,-4) = 



V4-CM 2 -C;V^ 2 +4 ^ figuration. ForC. 



where in addition to the light-induced dipole force Fj from 
Eq.(f3]), we include viscous friction with coefficient fi (related 
to the single-particle friction coefficient ha by ^ = VI^a)- For 
plastic beads immersed in water, /i follows from the Stokes 
law; for atoms in vacuum, it can represent some laser cool- 
ing mechanism. This equation is nonlinear, as its solution in- 
volves integrating (fl]i to obtain the electric field for the force. 
We proceed by linearizing Eq. <[8J around an equilibrium con- 
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<C A we have 



2(1 + Cf) 



Thus, at steady state, \A±\ \A^\ = \C\\ —,...,= 

\C N \, and \Bi\ =,...,= \B N \ = \Di\ =, ... ,= \D N \, i.e., 
the pump lasers fill the structure unattanuated. 

Now consider the steady state of N > 1 identical, purely 
dispersive trapped clouds, with d =,...,= £v = £ < 2/ A. 
Since at every cloud \Cj/Bj \ — \Bj/Cj\ = A, the phase slips 
are all equal: \i =,•••;= Xn — X- Thus the equilibrium 



N 
1=1 



where the matrix D is defined by 



D 



_d_ 

dxi 



F j (S n = 0,n=l..N). 



(9) 



(10) 
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configuration is an equidistant lattice, Xj — ^ — a- 1 -r 
(j — l)d. The lattice constant d is clearly independent of N, 
and a decreasing function of the phase shift x - see Fig.[T] and 
the introduction of J3l -, explicitly 
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(n-xiCA)). 



(6) 



Stability analysis requires finding the eigenvectors of D and 
determining their dynamics. Details of this calculation are in- 
volved and will be published elsewhere, we outline the proce- 
dure below. The key tool is the transfer matrix (TM) method, 
as used in [4]. The TM of the whole optical lattice is a prod- 
uct of the TM's of a single block of the lattice, which consist 
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FIG. 3: (color online) Real (thick) and imaginary (thin) part of the 
first few eigenvalues z\ (continuous purple), Z2 (slashed blue), Z3 
(dotted green), for a lattice of N = 100 (a) and N = 10 (b) clouds 
of polarizability £ = 0.1 each. 



of the BS transformation followed by free propagation over 
length d. Since losses are neglected, the two eigenvalues of 
the TM of a single block are e ±lB with £ C. The parame- 
ter 0, related to the quasimomentum, is given by the solution 
of cos 6 = cos kd — £ sin kd |4]. We solve this equation ex- 
plicitly and obtain the surprisingly simple result 



sin6 = C4/2; vr/2 < < tt. 



(11) 



We next apply the TM method to a perturbed OL where the 
Tth cloud is displaced by an infinitesimal amount. The cal- 
culations lead to explicit formulas for the matrix D which we 
omit here for the sake of brevity. Two important properties 
of D must be mentioned. First, Dji depends only on / — j: 
D is Toeplitz matrix. In particular, for Djj < all clouds 
are trapped in identical wells. Second, D is not symmetric. 
This shows that Fj is not a conservative force: if it were, 
Fj = —d/dxj V(x\...xn) would imply that D is a Hessian 
matrix, symmetric by Young's theorem. Note that reflection 
symmetry of the system is broken by the pump asymmetry. 

The eigenvalue problem of a nonsymmetric real matrix is in 
general not trivial. We have found, however, that a generalized 
Fourier transformation with complex wavenumbers diagonal- 
izes D exactly. The analytical formulas for the eigenvectors 
Vb and eigenvalues Zb of D, with b = 0, . . . , N — 1, read 



[v b ]j = (Ve 2 ™ b y/ N 

z b = p Vv cose 



4 K/V sin 2 



(12) 



, (13) 



where j3 — 8k(I /c is related to the oscillation frequency 
ojq of a single cloud in a symmetric (incident laser intensities 
Iq = Ix = cq I -Eq 1 2 c/2) trap by muj 2 = (3. Due to the pump 



asymmetry the eigenmodes ( fl~2T > of the lattice are complex, ex- 
cept for 6 = 0, which is a distorted center-of-mass mode and, 
if TV is even, b — N/2, the density wave of highest wavenum- 
ber possible (jr/d). These two modes are always stable, as 
z n/2 < zq < 0. Since D is real, all other modes form con- 



jugate pairs: Zb 



'N-b 



and Vfa = v 



N-b- 



We briefly discuss 



the meaning of these eigenmodes below. 

Consider a pair of complex eigenvalues Zb — z* N _ h with 
< b < N/2, and the corresponding eigenvectors = v£. 
Both Re(vb) and Im(vb) describe density waves of wave- 
length Nd/b, modulated so their amplitude increases towards 
the stronger pump. Now time evolution by (0 does not lead 
out of the subspace of spanned by these modes: for any 
superposition £ = pRe(vf,) + qlm(vb) with K, Eq.© 

is equivalent to a single complex homogeneous second-order 
linear differential equation, whose general solution is 



Here c± = p± + iq± are arbitrary constants, and 



(k± + iuo±) 



-fx ± y/fJ, 2 + Amz\ 
2m 



(14) 



(15) 



with k_ < k + to fix notation. This corresponds to two su- 
perimposed density waves of wavelength Nd/b, one coprop- 
agating with the stronger beam (u>_ < 0), and one counter- 
propagating (ci> + > 0). Their phase velocities are given by 
Nd \u>±\ /(27T&). The copropagating wave is exponentially 
damped with constant k_ < 0, but the counterpropagating 
wave can be either damped or amplified. Thus, this pair of 
modes is stable if k + < 0, which corresponds to 



m(Imzf,) 2 < — /j^ReZb- 



(16) 



For symmetric pumping ^1 = 0, the matrix D is symmet- 
ric, its eigenmodes ( Tl2l are the Fourier components, and the 
eigenvalues ( Tf3l are all real and negative, thus the lattice is 
stable. Almost all modes have the same frequency as a sin- 
gle trapped cloud, z\ = z-x =,...,= zjy = —(3, except the 
center-of-mass mode, with zq = —j3J (1 + N 2 ( 2 ), which be- 
comes soft if N — > 00. 

With the introduction of a pump asymmetry A > 0, the 
eigenmodes and the eigenvalues acquire imaginary parts, and 
as A is increased, the real parts of the eigenvalues turn positive 
one by one. The first few eigenvalues are shown as functions 
of A for two examples in Fig. [3] In the "strong collective 
coupling", N( ^> 1 limit (Fig. [3] a), we observe clearly sepa- 
rated resonances. In this limit, whenever tt — < ir/N, we 
have \fv w 1, and the denominator of dT3l is approximately 
1 — sin 2 0/ sin (nb/N), which, with Eq.dTTb. places the res- 
onance for mode b at A ~ A b = 2bw/(N£). We remark that 
A = 2ir/(N() fits the boundaries between the shaded green 
areas of Fig. |2] almost perfectly for A < 1. Outside of the 
strong collective coupling regime (Fig. |3]b), the resonances 
are not well resolved. It may even happen (as in the plot- 
ted example) that mode 6 = 2 becomes absolutely unstable 
(Rez2 > 0) at lower A than mode 6=1. This causes the 
"shoulder" in the N = 10 instability limit on Fig. [2] At the 
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FIG. 4: (color online) Time dependence of position distortions £ 
(in color coding, in units of 10 _3 A) in an asymmetrically pumped 
overdamped optical lattice of N = 100 clouds with polarizability 
f = 0.1, after excitation of mode Re(vi) at t = with amplitude 
10~ 3 A. The continuous grey contour line is £ = 0. In a), the sys- 
tem is subcritical: A = 0.632, and = —0.55 — 6.88i. The 
excitation results in a density wave propagating towards the stronger 
beam, and dying out. In b), at supercritical asymmetry .4 = 0.655 
the eigenvalue is zi//3 = 1 .48 — 5.94i. The density wave is now am- 
plified, and at t ~ 0.5/iAc/7o ~ 2.5/i//3 we leave the linear regime. 
Then a local drop in the lattice constant develops at x ~ 30A, which 
will result in two clouds coalescing, and eventually all particles will 
be pushed away by the stronger beam (not shown in figure). 



critical asymmetry A = 2/£ = 20, we have O — ir/2 and all 
eigenvalues are 0; for A > 20 all modes are unstable. 

A few remarks about the nature of these eigenmodes and 
the instability are in order. Two timescales govern the dy- 
namics of the OL: t q = ym/ |Rez{,| of the oscillations 
and Td = m/fj, of damping. For weak damping r Q <C Td, 



modes with nonzero Im zi, are potentially unstable, but damp- 
ing can restore their stability, cf. Eq.([T6il. At the other ex- 
treme, in the overdamped limit Td <C r D , the dynamics is ef- 
fectively first-order, and the copropagating mode disappears 
(is "damped out"), for the counterpropagating mode we have 
u) + = — Imzfc//i, and k + = Kezi,//!. Even with arbitrarily 
strong damping, the OL becomes unstable if Rc Zb > 0, as the 
rhs of dToT ) is negative. This "absolute instability" is used to 
define the shaded areas of Fig. [2] We illustrate the dynamics 
close to the absolute instability limit in Fig. [4] showing the 
results of numerical integration of Eq. tJH]) in the overdamped 
regime near this limit. 

Dynamical instabilities resulting from asymmetric pump- 
ing have been observed in a far-detuned OL where atom- 
light interaction was amplified by a ring cavity lfl3ll . In free 
space near-resonant light has to be used (detunings of a few 
tens of atomic linewidths seem realistic), and thus the influ- 
ence of spontaneous photons poses serious experimental lim- 
itations. We checked via simulation that the dissipative scat- 
tering force induces quantitative, but no qualitative changes as 
long as |Im£| < |Re£| /100. However, spontaneous emission 
also heats the clouds, putting an upper limit on the timescale 
accessible by an experiment, and complicating the very cre- 
ation of the OL. One possible way to circumvent the latter 
problem could be creating the OL at larger detuning, where 
spontaneous heating is negligible, and then continuously de- 
creasing the detuning of the trap beams down to the desired 
value. As for the timescale of an experiment, we estimate 
that, e.g., for a cold gas of Rb atoms in a dipole trap detuned 
by A = -IO7. . . - 2O7, forming N = 100 disk-shaped 
clouds, at pump power ratio V = 10, the destabilization rate 
k + can exceed the heating rate by orders of magnitude if the 
surface density of the clouds is t] > 1/(2A 2 ). 
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